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METHOD AND DEVICE FOR ASSIGNING WEIGHTING COEFFICIENTS FOR 
PERFORMING ATTITUDE CALCULATIONS WITH A STAR SENSOR 



General technical domain 

This invention relates to the domain of methods for 
determining the orientation of objects in space. 

More precisely, it relates to a method of assigning 
weighting coefficients to measurements of a succession of 
stars acquired by a star sensor in order to determine a 
spatial orientation. 

State of the art 

An increasing number of artificial satellites are 
orbiting around the Earth or are launched into space. 
Obviously, it is important to know the position of satellites 
in space, but it is also essential to know their attitude, in 
other words their orientation with respect to an inertial 
coordinate system in the celestial dome. 

A star sensor is an instrument that supplies its own 
attitude measurement calculated from direction measurements of 
stars that it carries out, or that directly transmits the list 
of measurements of elementary directions to a client device, 
for example such as the central software on a satellite on 
which the sensor is located. 

Figure 1 shows a satellite 1, for example in an orbit 2 
around a planet, for example the Earth, at an altitude 13 
above the planet surface. The satellite 1 comprises a star 
sensor 4 on one of its walls facing the stars in the celestial 
dome, and with a view field diagrammatically shown by the 
trunk of the. cone 5. The circle 8 represents the part of the 
celestial dome observable in the view field 5. The stars 10 
and 11 diagrammatically represent some of the stars that can 
be observed by the sensor 4 . It can be understood that a 
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plurality or a large number of stars can be observed at the 
same time within the circle 8. The stars 9 and 12 represent 
stars in the celestial dome outside the view field of the 
sensor . 

5 The satellite 1 displaces in the orbit 2 along the arrow 

6. Consequently , it can be understood that the circle 8 also 
displaces on the celestial dome along the arrow 7. Thus, the 
star 9 used to be in the circle 8, and the star 12 will be in 
the circle 8 in the future of the satellite displacement. 

10 Most frequently, the attitude of the sensor 4 or the 

satellite 1 is calculated using only some of the stars present 
in the view field 8 of the star sensor 4. For example, about 
ten stars may be chosen to make an attitude calculation. The 
direction of the star 10 or 11 with respect to a coordinate 

15 system 15 related to the sensor 4 is measured in the attitude 
calculation. Since the direction of the star 10 or 11 with 
respect to a known inertial coordinate system 14 is also 
known, the direction of the coordinate system 15 (and 
therefore of the satellite 1 or the sensor 4) with respect to 

20 the coordinate system 14 of the celestial dome may be deduced. 

The precision of the attitude estimate depends on the 
choice of stars to be measured for the direction calculation. 

The choice of stars is also important when star direction 
measurements are transmitted directly and are used by the 

25 client device of the star sensor. 

In some cases, the stars may be associated with a weight 
instead of being selected or rejected, which weights the 
importance to be assigned to them in subsequent processing. 

In other cases, the star sensor itself makes a first 

30 selection of stars, makes measurements on these selected stars 
and may send them accompanied by a weight or transformed into 
a global attitude estimate. 
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The measurements may also originate from different star 
sensors onboard the same satellite. 

In all cases, a selection is a set of potential or 
previously made measurements of star directions, for which the 
5 validity date is identical. As we have already said, weights 
may be associated with each measurement. Over time, a star 
sensor 4 generates a succession of selections, measurements 
being made on all or some of the selected stars. 

Traditionally, stars chosen as being the best stars will 

10 be selected for as long as possible. In the early dates of 
stars sighting, it was difficult enough simply to lock onto a 
star. Therefore, once it was found it was not released. 

One of the criteria for choosing stars in a selection is 
particularly the magnitude of each star that can be observed 

15 in the view field. Thus, the brightest stars should be 
chosen, namely stars with the lowest magnitude. One of the 
other selection criteria is the distance of the star from the 
optical axis of the view field. More precisely, Figure 2 
shows the circle 8 moving along the arrow 7 across the 

20 celestial dome. Choosing stars remote from the optical centre 
83 and located within a reference area 81 within the circle 8 
can improve the precision of the estimate depending on the 
optical axis 83 of the sensor. The axis 83 is the most likely 
to be affected by measurement errors. Stars located in an 

25 area reference 82 will be close to the optical axis 83 at one 
time or another and will provoke measurement errors. A star 
that is too close to the optical centre 83 will be rejected 
and another star further from the optical centre 83 will then 
be selected. 

30 At the moment, selection methods give excessive 

importance to these two criteria, with the result that stars 
are not changed frequently. 

However, the above methods have disadvantages. 
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Some stars, for which the direction is particularly badly 
estimated by sensor 4, disturb the global attitude estimate 
throughout the duration in which they are selected. This 
duration may be fairly long when the attitude of the star 
5 sensor 4 varies only slightly and when its view field is 
fairly wide. Figure 3 diagrammatically shows such a 

situation. Figure 3 shows the error on the attitude of the 
sensor 4 as a function of time. Changes 30 in the error level 
of the graph are due to star selection changes for 

10 determination of the attitude. It is seen that errors due to 
the choice of a selection are relatively long compared with 
the oscillations 31 on each plateau of the curve, which are 
due to observation errors of each star in the selection. It 
is observed that the errors 31 cancel out due to averaging 

15 over a relatively short time compared with the time during 
which each selection is observed. In other words, the noise 
due to oscillations 31 can easily be filtered by existing 
techniques for processing of data received from star sensors 
by client devices, since it is within the high frequencies of 

20 the frequency spectrum of the signal acquired by the sensor. 
The fact that the selection changes only infrequently results 
in low frequency noise that is difficult to filter. 

More generally, methods of assigning a weighting 
coefficient according to prior art are incapable of 

25 controlling the noise dispersion with time due to each 
selection . 

Furthermore, some phenomena (for example distortion) vary 
depending on the position of stars in each selection in the 
view field. 

30 More generally, methods for assigning weighting 

coefficients according to prior art are incapable of 
controlling these, space-time dispersion phenomena for stars to 
which high weights are assigned. 



Presentation of the invention 

The invention proposes to overcome these disadvantages. 

One of the purposes of the invention is to divulge a 
method for assignment of weighting or weighting coefficients 
for a calculation of a spatial orientation in order to control 
measurement errors of a star sensor. 

Another purpose of the invention is to vary the 
dispersion of stars for which the direction measurements are 
used either by the star sensor itself or by the client device 
of the star sensor. 

Another purpose of the invention is to propose a method 
of enabling space-time dispersion of star selections. 

Finally, another purpose of the invention is to propose a 
method that is capable of taking account of previously made 
selections in order to spread low frequency noise related to 
star selections more or less throughout the spectrum. 

To achieve this, the invention proposes a method for 
assignment of weighting coefficients to measurements of a 
succession of stars acquired by a star sensor related to a 
client device to determine a spatial orientation, 
characterised in that higher or lower preference is given to 
refreshment of the positions of measurements with the highest 
weights and / or stars on which these measurements are made by 
the star sensor and / or its client device, so as to displace 
part of the power of the error associated with the set of star 
measurements within the frequency spectrum. 

This invention also advantageously includes the following 
characteristics, taken alone or in any technically possible 
combination : 

- in the calculation of the weights of measurements in a 
current selection, the reinforcement or attenuation takes 
place as a result of applying a distance weight associated 
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with each measurement in the current selection and 
characteristic of an average distance between firstly the said 
measurement and secondly the measurements in the previous 
selections and the other measurements in the current 
5 selection; 

the distance weight associated with the current 
selection measurement is calculated as a weighted average of 
the corresponding distances between firstly the said 
measurement, and secondly the previous selection measurements 
10 and the other measurements in the current selection 
respectively; 

- the weighting coefficient associated with the distance 
between a first measurement in the current selection and a 
second measurement in a previous selection or another 

15 measurement in the current selection includes a memory 
coefficient associated with the said second measurement, and / 
or the weight of the second measurement if it belongs to a 
previous selection or a temporary weight if it belongs to the 
current selection; 

20 - the distance calculation combines the angular distance 

between the two measurements, and an identity distance that 
depends on the difference in nature of the two stars for which 
the measurements are being made; 

- the memory coefficient of a measurement m± at a time t 
25 is defined using the following formula: 

Mem(m±/t) = \ixn^ [t ~ T<m i )] , where 

• T(m k ) is a validity date of a measurement m k 

• p and n are constants. 

a charge is assigned to each star for which a 
30 measurement is made, the charge summarising the weights 
assigned to the measurements made on the said star in the 
past, attenuated by the passage of time; 
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- the charge of the star e p is defined at an instant T by 
the following formula: 

Cha (e p , T) = ^[Ai/n^ x Mem^/T)] 

Einti )=e p 

where Mem(m±/T) is the memory coefficient of the 
5 measurement m± at time T, E(mi) is the star on which the 
measurement m ± is made, and A (mi) is the measurement weight m±; 

the charge associated with a star to which a 
measurement in the current selection is related is updated 
before it is used in the calculation of the weight associated 
10 with a measurement, using a coefficient that depends on the 
difference A between the current date and the last update date 
for this charge; 

- the coefficient may be a factor and is in the form n~ A , 
where n is a constant; 

15 - the coefficient may be additive and is in the form - pxA , 

where p is a constant; 

after calculating the weight associated with a 

measurement in the current selection, the charge associated 

with the star for which this measurement was made is updated; 
20 - the update is made by adding the weight associated with 

the measurement ; 

- a random function is used in the calculation of the 
weights; 

- the calculation of the distance weight is iterated with 
25 a temporary weight for measurements in the current selection, 

the distance weight being used to calculate a new weight 
itself used to calculate a new distance weight and so on, 
until convergence towards a final weight; 

- the digital values of the method are saved in memory 
30 and processing means of the sensor and / or the client device; 
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the renewal rate of stars with a large weight is 
increased by increasing the frequency of measurements of the 
star sensor and / or the client device; 

- the dispersion of the complete new selection is used 
5 directly in the weights, using processing means related to the 

sensor and / or client device; 

- processing means related to the sensor and / or the 
client device comprising a neurone structure are used to 
directly affect dispersion in the weights. 

10 The invention also relates to the device for implementing 

this method. 

Presentation of the figures 

Other characteristics, purposes and advantages of the 
15 invention will become clear from the following description 
that is purely illustrative and is in no way limitative, and 
that should be read with reference to the attached drawings in 
which: 

Figure 1, already commented upon, diagrammatically 
20 shows an artificial satellite in orbit about a planet, and 
comprising a star sensor; 

- Figure 2, already commented upon, diagrammatically 
shows a view field of a star sensor on the celestial dome; 

Figure 3, already commented upon, diagrammatically 
25 shows a graph of the error on the attitude as a function of 
time in methods for assignment of weighting coefficients 
according to the state of the art; 

- Figure 4 diagrammatically shows a graph of the error on 
the attitude as a function of time in a method for the 

30 assignment of weighting coefficients according to the 
invention . 
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Detailed description 

In some of its calculation phases, a star sensor 4 shown 
in Figure 1 chooses to select specific stars (10 and 11 on 
Figure 1) that it detects within its view field 8, either to 
5 use it itself in the calculation of the attitude, or to 
transmit the measurement of its direction directly to the 
client device, for example the satellite 1. 

The selection then involves the assignment of a weight to 
each star 10 and 11 that characterises the potential of the 
10 star with respect to the needs of the client device 1 of the 
star sensor 4 . 

According to other methods, the stars are used in the 
calculation of the attitude with weighting, or a weight is 
transmitted to the client device 1 with the measurement of the 
15 star direction. All stars may then be transmitted, and a star 
transmitted with a zero weight is equivalent to an unselected 
star . 

In any case, a weight is associated with a measurement of 
a star direction at one time or another. This weight 
20 characterises the potential of the measurement to satisfy the 
needs of the client device 1 of the star sensor 4. 

We will now describe a first possible method for 
assigning weighting coefficients to measurements of a 
25 succession of stars acquired with a star sensor in order to 
determine a spatial orientation. According to this first 
method, the weights are assigned to the measurements 
individually. 

According to a first method, the importance assigned to 
30 measurements is controlled as a function of whether or not 
they relate to new stars, which may be called the time 
dispersion of star measurements. It takes account of previous 
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selections. The dispersion effect of measurements is 

amplified. 

According to a first possible embodiment, the consecutive 
duration during which a particular star is selected is reduced 
5 during observation of the stars. The effects of such a 
situation are shown diagrammatically in Figure 4. Thus, the 
error on the attitude of the sensor 4 in Figure 1 as a 
function of time is plotted in Figure 4. As in Figure 3, the 
changes 30 in the error level on the graph are due to changes 

10 in star selections made to determine the attitude. It is 
found that the errors due to the choice of a selection are 
relatively shorter than the errors in Figure 3. They are of 
the same order of magnitude as the errors 31 due to 
observation errors of each star in the selection. In terms of 

15 the spectrum, a higher frequency of changing the star 
selection is equivalent to replacing low frequency noise 
(induced by stars made noisy by oscillations 31 and kept for a 
long time) by higher frequency noise (induced by a more 
frequent change of stars used in the attitude estimate) . More 

20 precisely, this is equivalent to spreading the power at very 
low frequency within a wider frequency range. This operation 
may be particularly interesting if the client device of the 
star sensor has higher frequency measurement capabilities 
(which is the case when it uses high performance gyroscopes). 

25 Noise can be easily filtered. 

Conversely, it can be understood that if high frequency 
filtering capabilities of the client device are weak, and if 
the high frequency noise is bothersome, it may be useful to 
give preference to low frequencies, by keeping each star for a 

30 longer period. Therefore, there will be a mode for using the 
method controlling a different dispersion of measurements in 
time. 
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According to another possible method according to the 
invention, the dispersion of star measurements in space will 
be controlled. 

According to a first possible embodiment, when stars are 
5 being observed, the selected stars are geometrically dispersed 
in the view field. This controls measurement distortion. 
This dispersion relates to a selection made at a given 
instant . 

It can also be understood that according to another 
10 embodiment, it is required to not disperse the measurements in 
space . 

Space - time dispersion of measurements will be combined 
to control the error due to each selection in the orientation 
calculation . 

15 Consequently, in calculating the weight associated with a 

measurement, the invention uses the average distance of the 
corresponding measurement with measurements in previous 
selections and the present selection . If it is required that 
the measurements to which a high weight is assigned should 

20 last for a long time, then an attempt will be made to minimise 
this average distance . On the other hand, an attempt will be 
made to maximise the distance if it is required that the 
measurements should be renewed frequently. 

25 Definitions 

A "distance weight" means an average distance of a 
measurement with measurements in previous and present 
selections. The distance weight is then used in the 

calculation of the weight, which takes account of other 

30 considerations (for example the position of the measurement on 
the view field or the magnitude of the star considered) . One 
possible expression of the distance weight is given in the 
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remainder of this description- The first step is to define 
some of its possible variables. 

Let {mi,m 2 , ... m N } be the series of measurements of previous 
and present star directions* Two distinct elements in this 
5 sequence may correspond to the same physical star measured at 
two different instants. 

The sub-series {mi,rri2, ... m p }, where P < N, is the series of 
stars in the present selection. 

8(mi,mj) characterises the difference between measurements 

10 mi and rrij and uses two terms with very different natures, for 
example by their product: 

• DA (mi, rrij) - this is the angular difference in the 
coordinate system of the star sensor 4 (or any other 
corresponding distance, for example the geometric distance on 

15 the detection matrix of the same sensor) . Characterising this 
angular difference provides a means of dispersing or grouping 
measurements that will be assigned a high weight on the sensor 
matrix. For example, this provides a means of compensating 
for some optical distortion errors. 

20 • DE (mi, m-j) - this is the identity difference, in other 

words the difference in characteristics of stars for which two 
measurements are made. For example, this difference may be 
equal to 0 if the two stars are identical, and otherwise to a 
positive value (two measurements corresponding to identical 

25 stars can give different values of the distance DA, due to 
movement of the star sensor along the direction 6 in Figure 
1) . Characterising this distance makes it possible to replace 
stars for which measurements are made more or less frequently. 
This also provides a means of giving preference to stars that 

30 have spectra or magnitudes that are different from (or similar 
to) previously selected stars. 

One possibility consists of using the identity difference 
to characterise the difference between the current measurement 
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and measurements made in previous selections, and using the 
angular difference to characterise differences between 
measurements in the current selection. 

The distance may possibly be a function of one or both of 
5 the previous distances. If the function is decreasing both as 
a function of DA and DE, a proximity weight will be calculated 
instead of a distance weight. This will then be used in the 
formula for calculating the weight instead of or in addition 
to the distance weight. 
10 For example, we will choose 8(m±,mj) = qx (e DA(m j ,m i } - 

1) xD£(mj,irii) , where q is a constant, DA(mi,rrij) is the angular 
distance between the two measurements, and DE(rrii,mj) is equal 
to 0 if the stars for which the two measurements were made are 
identical, and otherwise 1. Therefore, S(m± r nij) is not zero 
15 only if the stars are different. 

The previously calculated weight A (mi) , which in 
particular uses the distance weight, is associated with each 
measurement m± in previous selections. 

The weight of previous measurements used in the 
20 calculation of the distance weight corresponding to a 
measurement mj of the current selection is reduced, because 
time has erased the memory of this measurement. The result of 
this reduction is the remanent weight. 

The remanent weight may be calculated as being the 
25 product of the weight and- a memory coefficient Mem(m±/t) that 
characterises the stored memory of the measurement m ± at date 
t. The measurement memory coefficient m± for the measurement 
mj at date T(rrij) is written as Mem (m±/T (nij) ) . 

The memory coefficient Mem(m±/t) decreases when the 
30 measurement mi becomes more remote in time. This reflects the 
gradual erasure of star measurements from memory. 

For example, the following product will be used: 

Mem(m±/t) = pxrr Et ~ rr V 3 ; where 
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• T(m k ) is the validity date of measurement m k 

• p and EI are constants. 

The values of constants may depend on the difference 
S(mi,mj). Thus, by varying jj, the memory coefficient may 
5 change to the value 0 when [t-T (m±)] exceeds a certain value. 
This erases the memory of the measurement. The coefficient 
may also change to the value 0 when [t-T(m±)] is less than a 
certain value. This provides a means of forcing the selection 
of stars present in very recent selections. It also provides 
10 a means of excluding the distance weight of measurements in 
the current selection from the formula. 

Calculation of the distance weight 

The distance weight associated with the measurement nij is 
15 an average over the set of measurements in previous selections 
and measurements in the current selection (except m-j). 

The average may be weighted by remanent weights of 
measurements used in the average. 

For example, the average will be calculated in the 
20 following form: 



P(m j ) = 



[A( m i )xMem(m i I T{nij )) J° x8 (m . , m i ) * 



^[Aim^xMemim^Tim^f 

1=1 



\/<p 



where 



Q and cp are constants. If Q and cp are equal to 1, then 
it is a linear average of the distances. If <p is equal to 2, 
it is a root mean square of the same distances (in principle Q 
25 is equal to 1). If <p < 1, the separation between pairs of 
very close measurement is accentuated. 

If the measurement nij is included in the present 
selection, the weight A(irij) will not use the distance weight. 
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In practice, the weight of a measurement in the current 
selection may be calculated from an arbitrary distance weight. 
An iterative method may be used after this first calculation: 
weights associated with these measurements in the current 
selection are used in the calculation of distance weights for 
these measurements, and take account of the distance weight 
calculated in the previous iteration. All distance weights of 
measurements in the current selection are thus recalculated 
until they converge. 

Another possibility consists of not using the other 
measurements in the current selection in the calculation of 
the distance weights of the current selection. This is 
equivalent to choosing Mem(m±/t) = 0 for all values of i such 
that m± is a measurement in the current selection. 

The calculation of the distance weight can give a weight 
for each star measurement. 

For example, the weight may be calculated as 
A(m j )=e~ M x x F(m j )xP(m j ) , where M is the magnitude of the star for 
which the measurement rrij is made, ^(mj) is the angular 
distance between the measurement irij and the optical axis of 
the star sensor, and P(irij) is the distance weight associated 
with the measurement m j . 

Thus, the stars with the greatest weights are given 
preference in a selection. 

It can then be understood that the space - time 
dispersion of selections is controlled by varying the 
different distance weight variables, which results in a 
different final weight to the measurements. The selections 
may be change more or less quickly, and the selections are 
more or less dispersed among themselves on the matrix of the 
sensor, as a function of the parameters of the chosen 
variables. 
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We will now describe a variant of the method in which the 
distance weight of the measurements in previous selections is 
calculated recursively, to limit the amount of information to 
be memorised and to simplify the calculations. 
5 In this variant of the method, the identity distance is 

used to define proximity between measurements. 

Definitions 

Depending on the variant of the method, the information 
10 that summarises the past characterises the stars on which the 
measurements in previous selections were made, rather than the 
measurements themselves directly. This information can be 
summarised as a pair C p associated with each star e p , for which 
the measurement formed part of a previous selection defined 
15 by: 

C p = {Cha, date}, where 

date is the validity date of the charge Cha, which varies 
with time. The charge Cha is the accumulation of remanent 
weights associated with the set of measurements at the instant 
20 considered, that in previous selections applied to star e p . 
Therefore the charge associated with star e p at time T is 
theoretically defined as follows, using the same notations as 
above and if E(mi) is the star to which the measurement mi 
applied: 

25 Cha(e p ,T) = £[i4(/w,) xMemCm/T)] 

i=P+\ 
E(m i )=e p 

The term inside the sum can be multiplied by a function 
of the angular distance between the current measurement 
applied to star e p and the measurement mi, to take account of 
the spatial dispersion between successive selections (the 
30 target area in the sky possibly changing under the effect of 
the movement of the satellite) . For example, this function may 
be in the following form: 
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_ -rixd (m (e ) ,m ) 
e pi 

where m(e p ) is the measurement in the current selection 
applicable to star p, d is the distance between the two 
measurements in the argument and r\ is a coefficient. 
5 Instead of a charge, a memory associated with the star 

may be calculated equal to the average of the memory 
coefficients associated with the measurements that applied to 
the star considered weighted by their weights, for example: 

N 

YMirn^xMemim.lT)] 
SOU(e p ,T) = ^ - 

E(ntt)=e p 

10 This memory is a "normalised charge" that will be used in 

the calculation of the weight. 

Propagation of charges 

Charges are propagated as follows at the time that the 

15 measurements in the current selection are processed. 

Charges are refreshed and associated dates are updated. 
Charges associated with stars that form part of previous 
selections are diminished to translate erasure from memory. 
This update may possibly concern only stars in the current 

20 selection, but all information related to other stars must be 
kept. Charges associated with stars that do not form part of 
any previous selection remain zero. Dates associated with 
charges of stars in previous selections change to the date of 
the current selection. 

25 This reduction in the charge may be made if A is the time 

difference defined as follows: 

A = current selection date - date previously associated 
with the charge, 

by using one of the following two methods: 
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1. multiplication of a previous charge by a factor: 

IT A , where n is a constant 

2. addition of a factor to the previous charge: 

-pxa, where p is a constant. 
5 The charge thus obtained is forced to 0 if the result of 

the previous operation is negative. 

The charge is calculated without approximation when the 
previous method has been chosen and the memory coefficient is 
also in exponential form. In other cases, the calculation is 
10 approximate. 

The next step is to assign distance weights. The 
distance weight assigned to each measurement in the current 
selection is equal to the refreshed charge associated with the 
star for which the measurement considered is made. If no 
15 previous measurements were made for the star, its initial 
charge is zero and therefore its distance weight will also be 
zero . 

The next step is to calculate the weights. The distance 
weights calculated in the previous step are used to calculate 

20 weights associated with measurements in the current selection, 
for example using the formula A (itij) =q~ m x y/ (itij) * P (itij) . 

The charges are then updated. Charges of stars to which 
measurements in the current selection are made are increased 
by the value of the corresponding measurement weights. For 

25 each star e h that was not previously characterised by a pair C h 
= {Cha, date}, the date is fixed at the current selection 
date. For all stars concerned, the date may be assigned in 
this step rather than in the charge refreshment step. 

Finally, charges that are too small are eliminated. This 

30 step is optional. Charges below a minimum value may be 
eliminated, which is equivalent to forgetting that at least 
one selection was made on the star in the past. 
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A random function, for example a random Gaussian or 
uniform variable may be used for calculating the weights. 
This contributes to space - time dispersion of stars, in this 
case by chance. 

5 

Using the method 

The distance weight or. the charge pair is calculated 
assuming that information is memorised about measurements in 
previous selections, for example chosen from the following for 
10 each measurement: 

• The measurement date, 

• The weight finally assigned to the measurement, 

• The direction of the measurement, identified by its 
angle in a coordinate system 15 related to the sensor 4, or by 

15 vector in a coordinate system related to its matrix, 

• Some characteristics of the star used to make the 
measurement . 

If the recursive method is used, the pair C p = {Cha, date} 
is kept for each star for which a set of measurements was made 
20 for which the memory was not erased. 

Therefore, memory and processing means are associated 
with the sensor and / or its client device, namely the 
satellite . 

If the available memory limit is reached, information 
25 about measurements for which the memory is weakest can be 
deleted. The criterion can apply to: 

• the product of the memory associated with the 
measurement by the weight of the same measurement; 

• the charge, particularly if the recursive method is 

30 used. 

We will now describe a second possible method for 
assigning weighting coefficients to measurements of a 
succession of stars acquired by a star sensor in order to 
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determine a spatial orientation. According to this second 
method, the best selection of stars for determining the 
spatial orientation will be chosen directly. We will no 
longer consider stars individually to determine which have the 
5 greatest coefficients before grouping them in selections. The 
dispersion of a complete new selection will be used directly 
in the weights, using processing means connected to the 
renewal rate of stars with a large weight is increased by 
increasing the frequency of measurements of the star sensor 

10 and / or the client device; 

According to this method, the processing means connected 
to the sensor may comprise a neurone structure. The neurone 
structure can assign weights to stars to determine the 
orientation of the sensor or the satellite, for example after 

15 a learning process. The learning criterion on the neurone 
network can use distance coefficients, average distances, 
charges and / or memory defined above. It may also use the 
average of distance weights, possibly weighted by the weights 
of stars (or by temporary weights that do not take account of 

20 distance weights) for the stars in a selection. 

In both methods, the selected star dispersion phenomenon 
may be amplified by forcing the star sensor to operate at a 
frequency higher than is strictly necessary. Provided that 
the previous methods are used, the number of selections used 

25 will be increased, causing additional spreading of noise on 
the spectrum. The advantages thus obtained could compensate 
the increase in noise induced by shortening of the integration 
of light information by the star sensor, from the point of 
view of the client device in the star sensor. 

30 Star measurements can then be filtered and then sub- 

sampled before use, to return to the required frequency. In 
principle, filtering involves the calculation of the attitude 
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quaternion. The method is efficient even if the previous and 
proximity weights are not used in the criterion. 

The above developments are equally applicable firstly to 
a star set detected and selected regularly (for example once 
5 every second) , at the same time (as is the case for 
conventional operation of a star sensor with a CCD detector) , 
or secondly to a star set in which the stars are detected 
irregularly (use of new types of detectors) and in which 
selections are produced irregularly (as is the case for APS 

10 (Active Pixel Sensor) star sensors. 

The methods may be used at several levels, possibly 
simultaneously, for example the selection level of the stars 
on which the stellar sensor (s) will make their measurements; 
weighting of measurements made at the star sensor (s); 

15 weighting of measurements made at the measurements client 
device . 



